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ABSTRACT 

We present a study of the hot gas and stellar content of 5 optically-selected poor galaxy clus- 
ters, including a full accounting of the contribution from intracluster light (ICL) and a com- 
bined hot gas and hydrostatic X-ray mass analysis with XMM-Newton observations. We find 
weighted mean stellar (including ICL), gas and total baryon mass fractions within rsoo of 
0.026 ± 0.003, 0.070 ± 0.005 and 0.096 ± 0.006 , respectively, at a coiTesponding weighted 
mean M500 of (l.OS^Q'^g) x lO^"' M©. Even when accounting for the intracluster stars, 4 out 
of 5 clusters show evidence for a substantial baryon deficit within rsoo, with baryon fractions 
(/b) between 50±6 to 59±8 per cent of the Universal mean level (i.e. the remaining 

cluster having /b = 75±1 1 per cent. For the 3 clusters where we can trace the hot halo to rsoo 
we find no evidence for a steepening of the gas density profile in the outskirts with respect 
to a power law, as seen in more massive clusters. We find that in all cases, the X-ray mass 
measurements are larger than those originally published on the basis of the galaxy velocity 
dispersion (<j) and an assumed cr — Afsoo relation, by a factor of 1.7-5.7. Despite these in- 
creased masses, the stellar fractions (in the range 0.016-0.034, within rsoo) remain consistent 
with the trend with mass published by Gonzalez, Zaritsky, & Zabludoff (2007), from which 
our sample is drawn. 

Key words: galaxies: clusters: general - X-rays: galaxies clusters - galaxies: stellar content 
- cosmology: observations - galaxies: evolution. 



1 INTRODUCTION 

The hot gas and stellar content of collapsed structures are key el- 
ements in the baryon budget of the Universe. Conducting a full 
census of these baryons is particularly important given the likely 
effects of non-gravitational heating and radiative cooling (e.g. Voit 
& Bryan 2001) on the hot gas and star formation. Both these mech- 
anisms raise the entropy of the intracluster medium (ICM), either 
by heating it directly or else condensing out the lowest entropy gas 
from the hot phase, but each produces a different balance of hot 
and cold baryons. However, tracking all of the baryons is observa- 
tionally challenging, particularly in the lowest mass haloes, where 
feedback may be more influential. In such cases, hot gas may have 
been displaced to the group/cluster outskirts or even ejected from 
the (progenitor) halo altogether (e.g. McCarthy et al. 201 1). More- 
over, a substantial fraction of the cold baryons may be locked up in 
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intracluster stars, whose diffuse and low surface brightness emis- 
sion is very difficult to detect (e.g. Feldmeier et al. 2002; Zibetti 
et al. 2005; Gonzalez, Zabludoff, & Zaritsky 2005; Krick & Bern- 
stein 2007). 

The importance of the ICL in the baryon budget is the subject 
of current debate, with different methods giving a range of esti- 
mates for its contribution to the total cluster stellar light. For exam- 
ple, direct measurements by Krick & Bernstein (2007) suggest an 
ICL fraction of 6-22 per cent in the r band, within 25 per cent of 
the virial radius, while Gonzalez et al. (2005, 2007) estimate that 
the combined contribution of the brightest cluster galaxy and ICL 
amounts on average to 33 per cent of the total stellar light within 
r2oo, based on detailed I band imaging of 24 clusters. In contrast, 
McGee & Balogh (2010) find a ~50 per cent contribution, based on 
mapping of hostless type la supernovae. In one of the largest studies 
of its kind, Zibetti et al. (2005) favour a lower contribution, of ~11 
per cent within 500 kpc, based on stacking Sloan Digital Sky Sur- 
vey (SDSS) images of 683 clusters. However, the use of stacking 
to map the ICL will tend to bias the inferred measurements towards 
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Cluster name 




Dec." 


Redshift 


Hi Column* 


Obsid'= 


Observation Date 


tlvIOSl 


tMOS2 






(J2000) 


(J2000) 




(xlO^O cm-2) 






(ks) 


(ks) 


(ks) 


Abell2401 


21 58 22.5 


-20 06 15 


0.0571 


2.39 


0555220101 


2010 Oct 29 


54.6 


54.7 


41.8 


Abell 2955 


01 56 59.8 


-17 02 23 


0.0943 


1.77 


0555220201 


2008 Aug 02 


75.6 


75.5 


51.5 


Abell S0296 


02 46 51.5 


-42 22 41 


0.0696 


1.86 


0555220301 


2008 Dec 26 


50.7 


52.9 


34.1 



Table 1. Key properties of the 3 galaxy clusters with new XMM-Newton data. "Coordinates of the brightest cluster galaxy, on which the X-ray halo is centred. 
''Galactic Hi column, interpolated to the cluster centroid using the data of Kalberla et al. (2005). '^XMM-Newton observation identifier. The last 3 columns are 
the effective exposure times for the MOS and pn detectors after flare filtering. 



lower values, if uncollapsed systems are included in the sample. 
The relatively wide variation in ICL estimates may partly be due 
to differences in sample selection. For example, the Gonzalez et al. 
study preferrentially includes systems with a prominent brightest 
cluster galaxy (BCG). 

On the other hand, the contribution of hot gas is better deter- 
mined, and represents the largest component of the baryon budget 
in clusters. However, the best constraints are based on X-ray se- 
lected (and often archival) samples (e.g. Vikhlinin et al. 2006; Ar- 
naud, Pointecouteau, & Pratt 2007; Sun et al. 2009), which are bi- 
ased towards haloes with more centrally concentrated hot gas (e.g. 
Rykoff et al. 2008b; Dai et al. 2010; Balogh et al. 2011). Further- 
more, there is a long-standing tendency for X-ray studies to tar- 
get different clusters to those best studied at optical/(near)-infrared 
wavelengths. A full accounting of hot gas and stars (including ICL) 
in the same systems is therefore lacking, but is clearly essential to 
establish the overall balance and uncover any missing baryons. This 
is best achieved by securing X-ray observations of optically well- 
studied clusters, given the difficulty of mapping the ICL directly. 

The goal of this paper is, for the first time, to combine X-ray 
observations of the hot gas in low-mass clusters with a full galaxy 
+ ICL stellar mass analysis, using the former to measure both the 
gas and total mass. This approach calls for X-ray follow-up ob- 
servations of clusters for which the ICL analysis has already been 
completed, given the difficulty of directly mapping diffuse optical 
light. The low-mass scale is particularly interesting to study since it 
encompasses the regime where cluster properties become increas- 
ingly diverse and where there is a more equal balance between hot 
and cold baryons (e.g. Gonzalez et al. 2007). Our clusters^ were 
drawn from the optically-selected sample of Gonzalez et al. (2007, 
hereafter GZZ07), with the aim of targeting the lowest mass objects 
(cr < 500 km s~^) with the best quality optical data, for which the 
gas and total mass within rsoo could be measured by XMM-Newton 
in a reasonable exposure time (rsoo being the radius enclosing an 
overdensity of 500 with respect to the critical density). 

Throughout this paper we adopt the following cosmological 
parameters: Hq — 70 km s~^ Mpc~^ , Jim = 0.3 and JIa = 0.7, 
converting any quoted data from other studies to match this choice, 
where appropriate. Derived quantities scale with the Hubble con- 
stant as follows: total mass, M oc h^^; gas mass, Mgas oc hjQ^^\ 
gas fraction, /gas oc h^g^^; stellar mass, Afstar oc h^ff; stellar 
fraction, /star oc h^g. Our adopted value for the Universal baryon 
fraction is nthn/^rnhn = 0.169 ± 0.008 from WMAP7 (Jarosik 
et al. 2011). For the conversion of I band optical light to stellar 
mass, we assume a mass-to-light ratio of 2.65 (see Section 3.1). All 
errors are la, unless otherwise stated. 



^ We use the term 'cluster' to refer to all the galaxy groups/poor clusters in 
our sample 



2 X-RAY DATA ANALYSIS 

2.1 Data reduction and spectral profile deprojection 

The full sample of 5 clusters consists of three systems for which 
we acquired new XMM-Newton observations, plus a further two 
clusters drawn from the archive (see Section 2.3). Our three clusters 
were awarded time in XMM-Newton Cycle 7, and observed with the 
EPIC-pn operating in extended full frame mode and the EPIC-MOS 
in full frame mode, with the thin optical blocking filter. Observation 
dates and cleaned exposure times are given in Table 1 . A detailed 
summary of the XMM-Newton mission and instrumentation can be 
found in Jansen et al. (2001, and references therein). Reduction and 
analysis were performed using the XMM-Newton Science Analysis 
System (SAS V. 10.0.0). We used the XMM-ESAS analysis scheme, 
which is based on the background modelling techniques described 
in Snowden, Collier, & Kuntz (2004), which represents the best 
currently-available method for analysing XMM-Newton data. 

The EPIC-MOS data were processed using the EMCHAIN and 
MOS-FILTER tasks, and the EPIC-pn data using EPCHAIN and PN- 
FILTER. Bad pixels and columns were identified and removed, and 
the events lists filtered to include only those events with FLAG = 
and patterns 0-12 (for the MOS cameras) or 0-4 (for the pn). EPIC- 
MOS CCDs in anomalous states were excluded from further anal- 
ysis. These included Abell 2955 MOSl CCD 4 and Abell S0296 
MOSl CCDs 2 & 4 and M0S2 CCDs 4 & 5. Point sources were 
identified using the CHEESE-BANDS task, including all three de- 
tectors, and excluded from all further analysis. Sources associated 
with the centre of the extended emission were considered to be mis- 
identified cuspy cluster emission, rather than genuine point sources, 
and so were ignored. 

Since our goal was to determine the radial spectral profiles 
of the clusters, we extracted spectra in annular regions centred on 
the peak of the X-ray emission, which also coincides with the cen- 
troid of the brightest cluster galaxy (BCG). The version of ESAS 
included with SAS V. 10.0.0 is only able to produce particle back- 
ground spectra for the EPIC-MOS cameras, so our spectral anal- 
ysis was based on these cameras only. However, we note that the 
background is lower and better understood for the MOS than for 
the pn camera, which is advantageous for the analysis of faint, ex- 
tended emission. Moreover, we were able to trace the X-ray haloes 
of our targets to larger radii using ESAS, compared to a conven- 
tional MOS-l-pn analysis. 

The XMM-ESAS scheme attempts to subtract the majority of 
the particle component of the background using data taken with the 
telescope filter wheels closed, and model the remaining X-ray and 
particle components. Annuli were initially chosen to contain 2000- 
4000 net source counts, with a minimum bin width of 30", to ensure 
that the effects of blurring by the point spread function are negligi- 
ble. Where the source emission became too weak to meet this crite- 
rion, the outer region of the field of view was divided into 2-3 large 
annuli, chosen to contain a sufficient number of counts to provide 
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Figure 1. Combined 0.3-2.0 keV MOS/pn XMM-Newton images of tlie diffuse emission from the 3 newly observed clusters, with point-like sources removed 
and filled using the CIAO task DMFILTH. The images have been adaptively smoothed (with S/N = 20), particle background subtracted and exposure corrected. 
The position of the brightest cluster galaxy is marked with a '+'. Equivalent images of the 2 archival clusters can be found in Sivanandam et al. (2009). 



well-defined background-only spectra. Spectra were extracted us- 
ing the MOS-SPECTRA and MOS_BACK tasks, and grouped to con- 
tain at least 20 counts per bin. 

Spectra were fitted in XSPEC 12.6.0 following a procedure 
similar to that described in the XMM-ESAS cookbook". Energies 
outside the range 0.3-10.0 keV were ignored and all EPIC-MOS 
spectra were fitted simultaneously. An additional ROSAT All-Sky 
Survey spectrum was extracted in an annulus between 0.5-1° from 
the cluster centre, using the HEASARC X-ray Background Tool'^ 
and was also included to help constrain the soft X-ray background. 
The residual particle background in each MOS camera was mod- 
elled using a power law whose index was tied across all annuli, 
with normalizations tied between annuli in the ratio of the expected 
scaled average flux, estimated using the PROTON_SCALE task. Sim- 
ilarly, the instrumental Al Ka and Si Ka fluorescence lines were 
modelled using Gaussians whose widths and energies were tied 
across all annuli, but with normalizations free. 

The X-ray background was modelled as a combination of 
four components whose normalizations were tied between annuli, 
scaling to a normalization per square arcminute as determined by 
the PROTON-SCALE task. The cosmic hard X-ray background was 
modelled using a power law with index F = 1.46 and normal- 
ization = 8.88 X 10"'' photons keV"^cm"^s"^ at 1 keV. The 
local soft X-ray background was modelled by a 0.1 keV unab- 
sorbed APEC model representing the local hot bubble, an absorbed 
0.1 keV APEC model representing the cool galactic halo, and a 
0.25 keV APEC model representing the hotter component of the 
galactic halo or the intergalactic medium. The relative normaliza- 
tions of the APEC models were free to vary, and the temperature of 
the hottest component was also allowed to vary after the initial fit. 
The line of sight absorption was modelled using the WABS compo- 
nent, fixed to the galactic column density taken from the FTOOLS 
task nh (based on the data of Kalberla et al. 2005). The RASS spec- 
trum was fitted using only these X-ray background components. 

The contribution of the cluster was modelled using a depro- 
jected (using the PROJCT component) absorbed APEC model. Tem- 
perature, abundance and normalization were free to vary in each 



bin, and where the best fitting normalization approached zero (in 
the outer bins) it was fixed to zero. Abundances were measured rel- 
ative to the abundance ratios of Grevesse & Sauval (1998). Where 
values of abundance (or in some cases temperature) were poorly 
constrained, or where 'ringing' between bins was observed, pa- 
rameter values were tied between pairs of bins. The resulting best- 
fitting values of normalization and la eiTors were then used to cal- 
culate the gas density. Finally, we performed a cross-check of our 
background modelling in the outermost bin of A2955 (which is free 
of emission from the cluster halo), and found that the 0.3-lOkeV 
count rate predicted by the background component of the model 
was within 1.5 and 0.7 per cent of the measured rate for MOSl and 
M0S2, respectively. 



2.2 Cluster modelling 

In order to determine the total mass and gas fraction of the clusters, 
we fitted the deprojected gas temperature and density profiles with 
the phenomenological cluster model of Ascasibar & Diego (2008), 
following the method described in Sanderson & Ponman (2010). 
This is a physically-motivated, parsimonious yet flexible model, in 
which the ICM is modelled as a polytropic gas in hydrostatic equi- 
librium within a Hernquist (1990) gravitational potential, with a 
variable cool-core component. It has been shown to provide a good 
description of the temperature and density profiles of a wide va- 
riety of galaxy clusters, fully encompassing the mass range of the 
objects studied here. (Sanderson et al. 2009a; Sanderson & Pon- 
man 2010). We use mean molecular weight of 0.601 atomic mass 
units, based on a fully ionized plasma with a mean metallicity of 
0.5 Solar and using the Grevesse & Sauval (1998) XSPEC abun- 
dance table"*. The associated ratio of electron to hydrogen (strictly, 
the atomic mass unit) number density is 1.157, used in calculating 
the gas mass. A characteristic radius was assigned to each annulus 
using the emission-weighted approximation of McLaughlin (1999): 



0.5 



I ' out 



3/2 



2/3 



(1) 



^ ftp : / /xmm . esac . esa . int /pub/xmm-esas /xmm-esas . pdf 

^ http://heasarc.gsfc.nasa. gov/ cgi- bin /Tools/ ^ http://heasarc.gsfc.nasa. gov/ doc s / xanadu/xspec/ 

xraybg/xraybg . pi manual /XSabund . html 
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Figure 2. Gas temperature and density profiles for the clusters (error bars & hollow circles + eiTor bars, respectively) with the best-fit ADOS cluster model 
(solid line) and Icr error envelope. Also shown, in columns 4 and 5, are two additional clusters included to increase the sample size, as described in section 2.3. 
The innermost bin was excluded from the fit in all cases and is plotted with dashed error bars. Dotted curves indicate extrapolation of the model beyond the 
data and the blue dashed curve for A2401 is the (much worse) best-fit model including the penultimate bin and with the polytropic index free to vary. Note 
that for A2401 the outermost bin is excluded in both cases and so the corresponding density point has not been corrected for projection and is plotted as an 
upper limit (see text for details). 



The cluster model was jointly fitted to the deprojected gas den- 
sity and temperature profiles, excluding the central bin, since our 
focus is on global properties and we expect complex baryon physics 
on the scale of the BCG to invalidate the model assumptions here; 
the fitted data and model are plotted in Fig. 2. To improve the ra- 
dial coverage of the gas distribution, we also use the cluster model 
to determine a correction factor for estimating the density in the 
outermost bin. This point is often omitted on the basis that the emis- 
sion in the outer annulus has not been deprojected and so includes a 
contribution from gas distributed along the line of sight beyond the 
radius of the corresponding spherical shell. However, we can use 
the cluster model to estimate what fraction of the total emission in 
this annulus originates from the corresponding spherical shell. 

To do this, we consider a line of sight column located at the 
half-radius of the spherical shell, and intersecting its outer edge at 
two points. The correction factor is simply the ratio of the density 
squared integrated along this line within the shell to that evaluated 
between ± infinity - this weighting approximates the emissivity of 
bremsstrahlung as well as that of low-temperature line emission. 
This fraction is then multiplied by the density obtained by assum- 
ing all the flux in the outer annulus is associated with the volume 
of the corresponding spherical shell. With the outer density point 
thus determined, a new best-fit model is obtained and used to eval- 
uate a corresponding new estimate of the correction factor, until 
convergence is achieved. The final coiTection factors obtained are 
0.70 for A2955 and 0.73 for AS0296. This then allows the outer- 
most density point to be used in the analysis, albeit at the expense 
of introducing some model-dependency. 

For A2401, the two outermost annular bins are only partially 
covered by the XMM-Newton CCD chips, as a result of the off-axis 
position of the cluster halo (see left panel of Fig. 1 ), resulting in sys- 
tematically incomplete azimuthal coverage. For the outermost bin 
only 23 per cent of the annulus falls on the chip, while the penulti- 



mate bin has a coverage of 59 per cent. Recent work by Eckert et al. 
(2012) has found substantial azimuthal variation in ICM density 
in cluster outskirts, as a result of structural asymmetries and gas 
clumping, which could introduce significant systematic bias into 
the determination of the gas properties from incomplete outer an- 
nuli. Moreover, both of these bins include emission from beyond 
rsoo and even beyond r2oo for the outermost annulus, where the 
intracluster gas increasingly deviates from hydrostatic equilibrium 
(e.g Nagai, Vikhlinin, & Kravtsov 2007) and may also be subject to 
complex physics like clumping (e.g. see Urban et al. 201 1, and ref- 
erences therein), which would invalidate the model. We therefore 
exclude both the outer two bins from the fit and so are unable to de- 
termine a corrected density value for the outermost annulus, as this 
is model-dependent. However, since the uncorrected density repre- 
sents an upper limit, we plot it in Fig. 2 with an arrow extending 
from the la upper bound. 

The best-fit parameters for the final cluster models and their 
la errors (estimated from 500 bootstrap re-samples) are listed in 
Table 2, together with values of rsoo and the radius of the outer- 
most annular bin included in the fit. Following Sanderson & Pon- 
man (2010), we consider the residual fractional intrinsic scatter of 
the gas density data about the model (which measures the real vari- 
ance of the data about the model, beyond that expected from the 
measurement errors) as a useful indicator of the effectiveness of the 
model in describing the data (second last column of Table 2). This 
quantity is the mean value of [{yobs - Upredf - Verr] IVpred, 
with all negative values to zero, based on the observed values, pre- 
dicted values and measurement errors of y (see the appendix of 
Sanderson & Ponman 2010). For A2401 and A2955 this scatter is 
three and one per cent, respectively, in gas density- well below the 
levels of >5 per cent associated with merging/disrupted clusters for 
which the model performs less well (Sanderson & Ponman 2010). 
However, AS0296 has a 7 per cent residual intrinsic scatter in den- 
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Cluster 


To (keV) 


t 


a (kpc) 


a 




/ 


rsoo (kpc) 


rfi% (kpc) 


X^/dof 


p value 


/P 


fkT 


Abell 2401 


2.94 ±0.25 


0±0.11 


588 ± 135 


0.01 ±0.03 


0.73 


±0.18 


672 ± 20 


432 


13.3/4 


0.010 


0.03 


0.02 


Abell 2955 


2.9 ±0.30 


1 ± 0.00 


709 ±51 


0.01 ±0.01 


0.78 


±0.11 


656 ±32 


770 


13.0/4 


0.011 


0.01 


0.13 


Abell S0296 


3.91 ±0.38 


0±0.00 


750 ± 125 


0.01 ±0.02 


0.73 


±0.12 


772 ± 39 


997 


36.8/5 


0.000 


0.07 


0.06 


Abell 2984 


3.72 ± 1.06 


0.52 ±0.16 


870 ±382 


0.49 ± 0.48 


0.30 


±0.18 


738 ± 66 


861 


10.6/1 


0.001 


0.00 


0.23 


Abell S84 


7.69 ± 1.43 


0.35 ±0.07 


1337 ±337 


0.38 ±0.69 


0.18 


±0.08 


1054 ±85 


963 


52.4/3 


0.000 


0.06 


0.06 



Table 2. Best-fitting model parameters for each cluster and Icr errors estimated from 200 bootstrap resamples. rfjt is the radius of the outer annulus of the 
fitted data in the spectral profile. The last four columns are the I degrees of freedom with associated p value and the mean fractional intrinsic scatter in the 
gas density and temperature, respectively (as defined in Sanderson & Ponman 2010). The lower portion of the table refers to two additional clusters included 
to increase the sample size, as described in section 2.3. 



Cluster 


Af25Oo(lOl*M0) 


M5oo(lOi4Af0) 


^2500 (kpc) 


rsoo (kpc) 


/gas(<»'2500) 


/gas«?"500) 


a" 


Abell 2401 


0.34 ±0.05 


0.91 ±0.08 


282 ± 15 


672 ± 20 


0.047 ± 0.003 


0.073 ± 0.009 


-0.48 ±0.17 


Abell 2955 


0.28 ±0.06 


0.88 ±0.13 


260 ± 17 


656 ±32 


0.036 ± 0.002 


0.064 ± 0.005 


-0.14 ±0.01 


Abell S0296 


0.48 ±0.12 


1.4 ±0.21 


315±27 


772 ±39 


0.039 ± 0.003 


0.065 ± 0.008 


-0.29 ± 0.25 


Abell 2984 


0.37 ±0.07 


1.3 ±0.35 


286 ±20 


738 ±66 


0.087 ± 0.004 


0.096 ±0.014 


-0.32 ±0.10 


Abell S84 


1±0.18 


3.7 ±0.96 


399 ± 23 


1054 ±85 


0.080 ±0.015 


0.075 ±0.019 


-0.47 ±0.11 



Table 3. Derived quantities from the cluster models with Itr errors estimated from 200 bootstrap resamples. "Logarithmic gradient of the gas density profile at 
0.04r6oo (after Vikhlinin et al. 2007); a < —0.85 indicates a strong cool core (Sanderson et al. 2009a). The lower portion of the table refers to two additional 
clusters included to increase the sample size, as described in section 2.3. 



sity, which hints at a somewhat irregular gas distribution which the 
model does not capture well. The fractional intrinsic scatter for the 
gas temperature is also presented in Table 2, although this is a less 
reliable diagnostic of model effectiveness (Sanderson & Ponman 
2010). 

According to the best-fit models, none of these clusters hosts 
a sharply peaked cool core: the inner cuspiness as measured by the 
logarithmic slope of the gas density at 0.04r5oo (a, cf. Vikhlinin 
et al. 2007; Sanderson et al. 2009a) is flatter in all cases than the 
value of < —0.85 which typically characterizes strong cool core 
clusters (Sanderson et al. 2009a, see Table 3). On the other hand, it 
is noticeable from Fig. 2 that the (excluded) innermost temperature 
point (radius < 40-50 kpc) lies substantially below the extrapolated 
model prediction in the cases of A2955 and AS0296. Our Ascasibar 
& Diego cluster model is fully able to describe cool cores (cf. As- 
casibar & Diego 2008; Sanderson & Ponman 2010), but hydrostatic 
equilibrium is built into this model, and so deviations are expected 
on the scale of the central galaxy. In any case, the possible presence 
of small cool cores in these two clusters has no significant impact 
on the present study, which is concerned with the large scale prop- 
erties of these systems. 



2.3 Analysis of two extra poor clusters 

In order to increase our coverage of the baryon fraction in 
groups/poor clusters for which a full intracluster light analysis has 
been performed, we have added two extra systems to our sample 
from a companion paper (Abell 2984 & Abell S84; Gonzalez et al., 
in preparation). These objects were also part of the GZZ07 sample 
and their XMM-Newton data were analysed in Sivanandam et al. 
(2009). For consistency, we have re-analysed these poor clusters in 
the same way as for our other three systems (i.e. as described in 
Sections 2. 1 & 2.2), and their corresponding model parameters and 
derived quantities are listed in the lower portions of Tables 2 & 3 . 
The corresponding correction factors for determining the gas den- 
sity in their outermost bins are 0.83 & 0.82 for A2984 & AS84, 
respectively. We address the comparison between our modelling 



analysis of these systems and that of Sivanandam et al. (2009) in 
Sections 5.3 & 5.4. 



3 OPTICAL DATA 

Cluster membership was determined on the basis of galaxy veloci- 
ties, which were obtained from dedicated spectroscopic (and asso- 
ciated photometric) observations in the case of A2955 and AS0296, 
yielding a total of 22 and 34 cluster members, respectively (Zarit- 
sky, Gonzalez, & Zabludoff 2006) . These galaxies were originally 
selected for spectroscopic follow-up if they were both fainter than 
the BCG and located within a projected distance of 1.5 Mpc of it, 
and the resulting velocity distribution was subjected to a 3(t clip- 
ping algorithm, to select the cluster members (Zaritsky et al. 2006), 
using biweight estimators of location and scale for improved ro- 
bustness (Beers, Flynn, & Gebhardt 1990). 

The positions of galaxies are plotted in Figs. 3 & 4 with point 
size indicating their / band absolute magnitude and with cluster 
members colour-coded according to their relative line-of-sight ve- 
locity (hollow points have unknown velocity). In addition, con- 
firmed non-cluster members have been omitted- a total of 35 and 
1 1 galaxies for A2955 and AS0296, respectively. Overlaid are con- 
tours of surface number density which were computed by kernel 
smoothing the positions of all the galaxies in the field with a Gaus- 
sian of bandwidth 500 kpc at the redshift of the cluster, using the 
task KDE2d in the MASS package (Venables & Ripley 2002) in 
r'. For A2401, galaxy velocities were extracted from the NASA 
Extragalactic Database (NED), giving a total of 24 cluster mem- 
bers. These were cross-matched with the 2 Micron All Sky Sur- 
vey (2MASS) Extended Source Catalogue (XSC; Skrutskie et al. 
2006) to provide K band magnitudes and their positions are simi- 
larly plotted in colour in Fig. 5, together with those of other 2MASS 
XSC objects in the vicinity, for which the velocities are unknown 
(coloured in grey). 

In all cases the BCG lies very close to the peak of the number 

^ http://www.r-project.org 



© 2013 RAS, MNRAS 000, 1-17 



6 A. J. R. Sanderson et al. 




Figure 3. Positions of galaxies around Abell 2955 witli corresponding isodensity contours superposed, Gaussian smootlied on a scale of 500 kpc. The point 
sizes indicate the / band absolute magnitude (assuming all the galaxies are at the redshift of the cluster) and the filled points are the cluster members, coloured 
according to the velocity with respect to the cluster mean velocity, normalized by the velocity dispersion (see Table 4). The position of the brightest cluster 
galaxy is marked with a '+' and the dashed circle marks the field-of-view of the con'esponding XMM-Newton observation. 
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Table 4. Optically-derived data for the clusters. The last 4 columns are the optically-derived mass, velocity dispersion, number of galaxies with velocities and 
Bautz-Morgan type, repectively (from Zaiitsky et al. 2006, GZZ07). The lower portion of the table refers to two additional clusters included to increase the 
sample size, as described in section 2.3. 



density contours smoothed on a scale of 500 kpc. However, other 
localized galaxy concentrations are evident, particularly in the case 
of A2401, which suggest either the presence of substructures in 
the cluster or ongoing merging/disruption. Moreover, the density 
contours are based on all the galaxies detected along the line of 
sight, not just those with redshifts and that were judged to be cluster 
members. Therefore, some of the features seen could be the result 
of unrelated fore- or background structures. 

Fig. 6 shows the variation in line-of-sight velocity with dis- 
tance from the BCG for galaxies spectroscopically confirmed to be 
members of the cluster. The behaviour seen is broadly consistent 
with the expected distribution for a virialized population of tracer 
particles, with no major anomalies that might cast doubt on the use 



of (7 as a mass estimator. For AS0296, there are four galaxies in 
the top right corner of the plot which are likely to be interlopers 
lying outside the virialized volume of the cluster; indeed 3 of these 
are visible as a clump of red dots in the lower right corner of the 
galaxy distribution plot in Fig. 4, suggesting an infalling substruc- 
ture. However, their impact on the overall velocity dispersion of 
the cluster is minor and excluding their contribution would only 
further exacerbate the mass discrepancy compared to the larger X- 
ray determined value, discussed in Section 5.1. A2401 may also be 
affected by one or two interlopers, on the basis of Fig. 6, but here 
also any bias is likely to be small. 

We performed several tests of normality of the galaxy veloc- 
ity distribution, using the R package NORTEST (e.g. Thode 2002). 
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Figure 4. As for Fig. 3, but for Abell S0296. 




Figure 5. As for Fig. 3, but for Abell 2401 and using 2MASS K band magnitudes. In this case the cluster members (filled, coloured points) were selected 
using data from the NED and the ' X ' symbols mark the positions of 3 galaxies with velocities but no corresponding K magnitude. 



For all three clusters, none of the tests appropriate for these sam- 
ple sizes shows evidence for a non-Gaussian distribution, including 
the Anderson-Darling test, which has been shown to be the most 
sensitive to optical cluster substructure (Hou et al. 2009). 

A key element in determining the stellar mass in clusters is the 
aperture within which the measurement is made, which depends on 
the total gravitating mass. In the original GZZ07 analysis. X-ray 
data were generally not available and so the total mass (within rsoo) 



was estimated from the velocity dispersion of the galaxies, using 
the relation between a and M500 based on the cluster sample of 
Vikhlinin et al. (2006), where a full hydrostatic mass analysis was 
performed. There is evidence that a can often underestimate the 
total mass in galaxy groups (Osmond & Ponman 2004; Connelly 
et al. 2012), and so a key aim of the present study is to determine 
the mass independently, using our X-ray data. Indeed, as shown by 
Balogh et al. (2008), the four lowest mass GZZ07 clusters, includ- 
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Figure 6. Line-of-sight velocity vs. distance from tlie briglitest cluster galaxy for the galaxies spectroscopically confirmed as cluster members, point size 
coded by absolute magnitude. For Abell 2401 , ' X ' symbols mark the positions of 3 galaxies with velocities but no corresponding K magnitude. 



ing A2955 from our study, have uncomfortably high baryon frac- 
tions (implying an underestimated total mass) on the basis of their 
velocity dispersions-derived masses. However, these systems had 
been excluded by GZZ07 on the grounds that their mass estimates 
required an extrapolation of the i\/5oo — o" relation and were hence 
unreliable (we emphasize that they were not excluded on the basis 
of /b)- We address the differences between the X-ray and optically 
derived mass estimates in Section 5.1. 

3.1 Intracluster light and stellar mass 

GZZ07 made drift-scan observations to measure the optical light 
profile of the clusters, with accurately measured background and 
flat-fields, in order to map the diffuse ICL directly and with good 
precision in each case. The cluster radial light profile was modelled 
as a combination of two de Vaucouleurs profiles, to capture the 
separate contributions from the BCG and ICL. The total optical lu- 
minosity of the cluster is calculated by integrating the model within 
a given aperture and then adding in the luminosities of the individ- 
ual non-BCG cluster galaxies in the same aperture (corrected for 
contamination). Since our X-ray derived masses (Table 3) are sub- 
stantially larger than the masses derived from velocity dispersions 
by GZZ07, as we will discuss in Section 5.1, the apertures corre- 
sponding to r-2500 and rsoo, within which the optical luminosity are 
integrated are correspondingly larger. We therefore re-measured the 
cluster optical luminosity (including ICL) within these new radii. 

These aperture measurements represent the projected lumi- 
nosity within a cylinder along the line of sight through the clus- 
ter. We have therefore corrected them for projection, to calculate 
the 3 dimensional luminosity within a sphere, assuming a Navarro, 
Frenk, & White (1995) profile with a concentration parameter of 
2.9 (based on the stacked profile of Lin, Mohr, & Stanford 2004). 
This correction only applies to the galaxy distribution and not the 
BCG-l-ICL light, which is (approximately) completely contained 
within the smallest aperture considered here (r25oo)- The ratios of 
the deprojected to projected luminosities are in the range 0.74—0.91 
within r25oo and 0.78-0.86 within rsoo- 

Our adopted average / band mass to light ratio (M/L) is 2.65, 
which is lower than that used in GZZ07. This new value incorpo- 
rates a 15 per cent contribution from dark matter as well as a correc- 
tion to the previous GZZ07 M/L calculation which is described in 
more detail in a forthcoming paper (Gonzalez et al, in preparation). 
The corresponding stellar masses and mass fractions are summa- 
rized in Table 4, together with other optical data for the 5 clusters. 



In Section 5.5 we address the impact on our results of systematic 
errors in M/L. 



4 RESULTS 

4.1 Radial trends in stellar and gas fraction 

The radial distribution of baryonic components is plotted in Fig. 7, 
with the WMAP7 inferred Universal baryon fraction shown for 
comparison. The gas fraction from the cluster model is plotted as a 
solid curve with grey error envelope, with dots indicating extrapo- 
lation beyond the data. Aperture values of the stellar mass fraction 
(including the ICL) are plotted at r25oo and rsoo, with correspond- 
ing total baryon fraction (gas -I- stars) plotted at the same radii. It 
can be seen that for all 5 clusters the stellar fraction is higher within 
^2500 than rsoo, although this difference is not statistically signif- 
icant for A2401. This is also consistent with the observation that 
the ICL in these systems is more spatially concentrated than the 
galaxies themselves (GZZ07), as also found in the SDSS stacking 
analysis of Zibetti et al. (2005), for example. On the other hand, 
the gas fraction profiles tend to increase monotonically with radius, 
with the exception of AS84 and A2984, in accordance with the ten- 
dency for the intracluster medium to be more extended than the 
dark matter (e.g. David, Jones, & Forman 1995; Sanderson et al. 
2003; Vikhlinin et al. 2006; Sun et al. 2009). Interestingly, the 
nearly flat /gas profiles and more centrally-concentrated /star for 
AS84 and A2984 imply a decline in total baryon content with ra- 
dius. However, this difference is not significant for AS84 and, for 
A2984, /b within rsoo is only 1.7cr lower than /b within r25oo- 

The mass fractions for the baryon components averaged over 
all 5 poor clusters are listed in Table 5, measured within r2soo and 
rsoo and weighted by the inverse squared measurement errors. The 
calculation was performed in log space since the baryon fraction is 
right skewed, as a consequence of being strictly bounded by zero 
below. This is confirmed by the fact that a log-normal distribution 
provides a better description of the data than a normal distribution 
(as determined using the FITDISTR function in the MASS package 
in R; Venables &. Ripley 2002). It can be seen that the combined 
baryon fraction within the two radii does not significantly differ, 
reaching a value of 0.096±0.006 within rsoo (at the corresponding 
mean Msoo of (l.OStfS ig) X 1O"M0), and falls some way short 
(reaching only 57 per cent) of the Universal mean level of 0.169 ± 
0.008 from WMAP7 (Jarosik et al. 2011). On the other hand the 
balance of hot and cold baryons does vary radially, with the ratio 
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Figure 7. Gas mass fraction profiles plotted out to r200 for each cluster, with the stellar and total baryon fraction (including ICL) evaluated at rsoo and r2500 
(offset slightly in radius to avoid over plotting). The grey region is the Icr confidence envelope on the gas fraction profile from the cluster model and a dotted 
curve denotes extrapolation of the model beyond the radius of the outermost bin (see Table 2). 
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Table S. Weighted mean mass fractions for the baryon components and total 
mass, measured within rsoo and r2500. The calculation was performed in 
log space, yielding an asymmetric standard error on the unlogged values. 



of stellar to total baryon mass declining from 0.44 to 0.27 between 
r2500 and rsoo- 

Looking at Fig. 7, if we consider only the systems where the 
fitted data extend beyond rsoo (i-e- where no extrapolation is re- 
quired), it is clear that a substantial baryon deficit is present in 
two of three cases, and is a factor of ~2 for AS0296 (with a best- 
fitting value of 50 ± 6 per cent of the Universal mean); for A2984, 
/b = 0.127 ± 0.019, representing 75 ± 11 per cent of the Uni- 
versal mean /b. We demonstrate in Section 5.6 that, for 4/5 cases, 
the baryon fractions cannot be reconciled with the Universal mean 
even when conservatively accounting for systematic errors. 



4.2 Variation in stellar and gas fraction 

With our new X-ray derived total masses and the associated over- 
density radii, we are able to address one of the key results of 
GZZ07, namely the stellar mass fraction in poor clusters. The left 
panel of Fig. 8 shows /star vs. i\/5oo for these 5 clusters. Also plot- 
ted, as dotted error bars, are the positions of 4 of the clusters from 
the original GZZ07 analysis (the A2955 /star measurement was ex- 
cluded from that analysis, since its mass estimate was deemed unre- 



liable, being based on an extrapolation of the Afsoo - o" relation). It 
is clear that the new positions for these points are substantially dif- 
ferent, primarily due to the X-ray masses being considerably larger 
than those derived from the velocity dispersion, which we discuss 
further in Section 5.1. 

The larger masses result in a larger rsoo aperture within which 
/star is measured and, since the BCG and ICL stellar contribution 
falls off rapidly with radius (Section 4. 1), the result is a lower stellar 
fraction. However, it is interesting that the change in both quanti- 
ties acts to move the points mostly along the regression line, thus 
largely preserving the trend found by GZZ07. While this is not al- 
together surprising given the covariance between /star and Msoo 
(since /star oc A/^jq), the shift is nevertheless along a flatter tra- 
jectory than implied by the relation between these quantities. For 
reference, the dashed line in the left panel of Fig. 8 shows the best 
fitting orthogonal regression from an analysis of new observations 
of the (optically selected) GZZ07 sample (Gonzalez et al, in prepa- 
ration), which incorporates the results from this work. Despite the 
substantial change in total mass and a reduction in the measurement 
errors, the five clusters remain fairly consistent with the trend and 
scatter in the /star - Afsoo relation from GZZ07. 

The corresponding relationship between gas fraction and mass 
is plotted in the right panel of Fig. 8. For comparison, data from the 
X-ray selected cluster samples of Vikhlinin et al. (2006) and Ar- 
naud et al. (2007) as well as the X-ray selected groups sample of 
Sun et al. (2009) are also shown, with error bars omitted for clar- 
ity. The solid line is the best-fitting BCES orthogonal regression 
(Akritas & Bershady 1996) in log-log space, which excludes our 
5 clusters and thus represents a purely X-ray selected sample. Of 
the 5 clusters, only A2984 is located near the trend, with the re- 
maining 4 all lying near the lower edge of the scatter about the 
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Figure 8. Left: Stellar mass fraction (including intracluster stars) vs. cluster total mass within rsoo, updating the results from GZZ07 (dotted error bars, which 
were not corrected for projection) with X-ray derived data from this work (solid error bars). The dashed line is the best fitting orthogonal regression from an 
analysis of new observations of the GZZ07 sample (Gonzalez et al, in preparation). This sample includes the five clusters from this work, but mainly comprises 
more massive clusters. Note that A2955 was excluded from the original plot data of GZZ07, but would have been located off the left hand edge of the plot, and 
that an / band mass-to-light ratio of 2.65 is used to calculate stellai' mass (see Section 3.1 further details). Right: gas fraction vs. mass for the same 5 clusters, 
together with X-ray selected data from the literature (Vikhlinin et al. 2006; Arnaud et al. 2007; Sun et al. 2009). The solid line is the best-fitting BCES power 
law regression to the X-ray selected literature data only (i.e. excluding our 5 clusters), with logarithmic slope of 0.186 it 0.083. 



fit, in contrast to the ./star-Afsoo relation. This is not necessarily 
surprising, since these were not X-ray selected clusters and there- 
fore might be expected to have lower /gas (and hence Lx), as seen 
in other optically-selected groups and clusters (e.g. Bower et al. 
1997; Rasmussen et al. 2006; Rykoff et al. 2008a; Hicks et al. 2008; 
Dai et al. 2010; Balogh et al. 201 1). Indeed, Rykoff et al. (2008b) 
demonstrate that differences in the X-ray luminosity-total mass re- 
lation between optically and X-ray selected cluster samples can be 
at least partly explained by the combination of Malmquist bias and 
deviations from hydrostatic equilibrium affecting the latter. 



4.3 Variation in total baryon content 

The sums of the stellar (including ICL) and gas fractions within 
rsoo for the five clusters are plotted in Fig. 9, to show the varia- 
tion in total baryon content compared to the Universal mean frac- 
tion. Fractional errors on the total baryon fraction were calculated 
by adding in quadrature the fractional errors on the gas and stellar 
fractions. Also indicated in Fig. 9 is the level of a 10 per cent de- 
pleted Universal /b (i.e. Q.9xQ,i,/Q.m hereafter referred to as the 
'depleted cosmic' /b), representing the typical value seen within 
rsoo in simulated clusters (e.g. Ettori et al. 2006). It can be seen 
that 4 out of 5 of our (optically-selected) clusters fall significantly 
(> 2a") below this line, implying a baryon deficit. With the excep- 
tion of A2984, this is mainly the result of a low gas fraction com- 
pared to the trend (right panel Fig. 8), although A2984 lies close 
to the trend in both gas and stellar fraction but still shows evidence 
of baryon shortfall. Furthermore, even when allowance is made for 
systematic uncertainties, these baryon fractions cannot be recon- 
ciled with the Universal mean (see Section 5.6). 



5 DISCUSSION 

An important aspect of our findings is that the X-ray determined 
masses of the 5 clusters are all substantially greater than implied by 
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Figure 9. Total baryon fraction (hot gas + stars) vs. total mass within rsoo 
for the 5 clusters. The horizontal dashed line shows the 'depleted cosmic' 
/b typically expected within rsoo from simulations without radiative cool- 
ing or feedback (e.g. Ettori et al. 2006). 



the velocity dispersion of their galaxy members. Therefore, in this 
section we first address this issue, before turning to a discussion of 
potential systematic effects relating to the gas and stellar mass mea- 
surements and finally considering the interpretation of our results 
and their associated implications. 



5.1 X-ray vs. optical mass estimates 

Estimating the gravitating mass of galaxy clusters is crucial in two 
ways in determining their baryon fraction, since this enters as the 
denominator in that quantity as well as determining the aperture 
within which the measurement is made. The original GZZ07 anal- 
ysis used the optical velocity dispersion of the galaxies to estimate 
M500 (via the M — Tx and a — Tx scaling relations) and this quan- 
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Figure 10. A compaiison of M500 as measured from X-ray data in this 
work witli the optically-derived estimates in GZZ07, which were based on 
the velocity dispersion of galaxies in the cluster and an assumed it- A/500 
relation. The dashed line shows the locus of equality. 



tity is plotted in Fig. 10 against the corresponding X-ray inferred 
value measured in this work, with the locus of equality indicated. 
In all cases the X-ray estimates are significantly higher than their 
optical counterparts, by a factor of between 1 .7 and 5.7 (mean 3.1). 

Part of the discrepancy can be attributed to the calibration of 
the M — a used by GZZ07, which was necessarily limited by the 
lack of clusters with suitably well measured masses and velocity 
dispersions. However, other potential causes of such a discrepancy 
are failures of one or both of the two key assumptions in both X- 
ray and optical mass estimates, namely equilibrium and spherical 
geometry. On the X-ray side, non-thermal pressure support and/or 
deviations from hydrostatic equilibrium (e.g. due to recent disrup- 
tion) are expected to cause X-ray mass wnrferestimates of ~5-20 
per cent (e.g. Kay et al. 2004; Rasia et al. 2006; Nagai et al. 2007; 
Piffaretti & Valdarnini 2008); accounting for this effect would only 
exacerbate the disagreement in this case. Nevertheless, deviations 
from hydrostatic equilibrium are likely to contribute towards the 
scatter in X-ray mass estimates, for which we make allowance in 
our analysis of the systematic error budget in Section 5.5. 

With the exception of A2984 (where the uncertainties are 
large), none of the clusters hosts a strong cool core (see Table 3), 
the absence of which is often associated with recent disruption 
(e.g. Sanderson, Ponman, & O' Sullivan 2006). On the other hand, 
there is no indication of any significant offset between the X-ray 
peak/centroid and the BCG, which suggests a lack of disturbance 
(Sanderson et al. 2009a). Moreover, since these are optically se- 
lected clusters, it is plausible that they may have a lower incidence 
of cool cores compared to X-ray selected samples, since the pres- 
ence of a cool core invariably boosts the X-ray flux in the region of 
highest surface brightness. 

Another possible explanation is that the line-of-sight veloc- 
ity dispersion is systematically underestimated for these clusters, 
which could occur if the cluster was substantially elongated in the 
plane of the sky. This hypothesis is difficult to test, given the limited 
number of member galaxies, but the morphology of the isodensity 



contours of all the surrounding galaxies (which includes fore- and 
background objects as well as cluster members) hints at such a pos- 
sible elongation in the case of A2955 (Fig. 3) and AS0296 (Fig. 4). 
No such elongation is evident for A2401 (Fig. 5), and here the opti- 
cal and X-ray masses are in better agreement (although still signifi- 
cantly discrepant). Interestingly, Connelly et al. (2012) have found 
that dynamical masses have a large scatter and, for less massive 
clusters in particular, are biased low compared to X-ray mass esti- 
mates. This suggests that, for poor clusters in general. X-ray mass 
estimates may be more reliable than those based on the velocity 
dispersion of the member galaxies, particularly when the number 
of galaxies is relatively small (~20-30; Table 4). Indeed, Biviano 
et al. (2006) find that the scatter in cr-based mass estimates roughly 
doubles when the number of galaxies decreases from ~400 to ~20 
(reaching --^0.6 for 7Vgai=10), based on an analysis of 62 simulated 
clusters. However, this alone is not sufficient to account for the dis- 
crepancy in this case, which amounts to a systematic bias of at least 
a factor of ~2. 



5.2 The impact of the total mass on /b 

The total gravitating mass determines the aperture for measuring 
the gas and stellar masses, via the definition of an overdensity. A, 
with respect to the redshift-dependent critical density of the Uni- 
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The stellar and gas fractions are inversely proportional to A/a, by 
definition, and in general both vary oppositely with radius (Fig. 7). 
As a result, /b within a given overdensity radius, ta, depends quite 
subtly on the corresponding A/a, via Equation 2. Thus a large 
change in A/500 does not necessarily change the baryon fraction 
within rsoo substantially. As a consequence, /b is quite robust to 
variations in A/500 resulting from modelling uncertainties, for ex- 
ample. 

The sensitivity of the baryon fraction measured within r5oo to 
the assumed A/500 can be judged by exploring the effect of arti- 
ficially altering the value of r5oo, as follows. The gas and stellar 
masses are measured within each nominal value of rsoo and the 
corresponding /b is calculated by adding them together and divid- 
ing by the Adsoo calculated from Equation 2. The stellar mass is 
calculated from fstar, which is fixed at the spot values for r25oo 
and r5oo for r ^ r25oo and r ^ 7-500, respectively, and linearly 
interpolated in log-log space in between. Fig. 1 1 shows curves of 
the resulting implied baryon fraction measured in this way within 
each nominal rsoo aperture. 

Under the hypothesis that the true /b is 15 per cent (roughly 
the 'depleted cosmic' /b level of 0.9 times the Universal mean), 
the corresponding rsoo required to achieve this result is plotted 
as a dashed vertical line in Fig. 1 1 . For comparison, the observed 
rsoo is shown as a dotted line and ratios of the implied A/500 and 
rsoo to their observed counterparts are shown in the strip titles for 
each panel in the plot. With the exception of A2984, the masses 
would need to have been overestimated by ~ 2-4 times to achieve 
/b = 0.15, equivalent to an overestimate in mean temperature of 
~ 1.3-2.7 times. Thus the observed deficiency of baryons can only 
be accounted for by a factor of ~ 2-4 systematic error in the deter- 
mination of the total mass. Such a large error greatly exceeds both 
the measurement errors and our best estimate of the true systematic 
uncertainty, as we demonstrate in Section 5.3. 
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Figure 11. Baryon fraction profiles implied by assuming tliat each radius corresponds to rsoo, with dashed vertical lines marking the radius required to 
achieve a nominal 'depleted cosmic' value of /t, = 0.15 within rsoo, as compared to the position of the observed value of r'soo (dotted lines). Ratios of the 
implied A/500 and rsoo for the case of /i, = 0.15 to the corresponding observed values are indicated in the strip text for each panel. See text for details. 



5.3 Cross-check of X-ray masses using the M ~ Tx relation 

To test the accuracy of the masses determined from the cluster 
modelhng, we also estimated M500 from the mass-temperature 
(M - Tx) relation of Vikhlinin et al. (2006), based on a global 
mean temperature (T). We estimated T using Pgas('') and T{r) 
from the best-fitting cluster model, evaluating a weighted (by 
Pgas's/r) mean in the radial range (0.15-1) rsoo, determined iter- 
atively (since rsoo depends on Msoo, which depends on T). The 
comparison between Afsoo calculated by these two methods is 
plotted in Fig. 12, which shows good agreement in all cases and 
no indication of any systematic bias. This is consistent with the 
good agreement in rsoo found by Sanderson & Ponman (2010) be- 
tween the ADOS cluster model and the Vikhlinin et al. M — Tx 
relation, using direct extraction and fitting of spectra in an itera- 
tively adjusted aperture (from Sanderson et al. 2009b). Although 
the M — Tx masses estimated in this way are not completely inde- 
pendent from the ADOS model masses, this comparison neverthe- 
less demonstrates consistency between both the two methods and 
the cluster samples, in respect of hydrostatic equilibrium. Since the 
Vikhlinin et al. clusters were carefully selected to have very regular 
X-ray morphology and with only weak signs, if any, of dynamical 
activity, it is therefore reasonable to infer that our cluster model- 
derived X-ray mass estimates are reliable and not substantially af- 
fected by any departures from hydrostatic equilibrium and/or spher- 
ical symmetry within the hot gas haloes of the five clusters. 

A further comparison of the total cluster mass is available for 
A29S4 & ASS4 from the independent analysis of the same XMM- 
Newton data by Sivanandam et al. (2009), also using the Vikhlinin 
et al. (2006) M — Tx relation, but measuring the mean temperature 
directly via spectral fitting of the raw data. The values of Afsoo ob- 
tained by Sivanandam et al. are (in units of 10^* M0) 0.96 ± 0.098 
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Figure 12. A comparison of Msoo determined from the cluster model with 
that estimated from the M - Tx relation of Vikhlinin et al. (2006). The 
dashed line shows the locus of equality. 
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Figure 13. A comparison of the best-fitting AD08 and /3-modeI gas density profiles, witli the measured data plotted as hollow points and (mostly very small) 
error bars. Within each panel the profiles extend out to r200- For A2401, the excluded outermost two points (see Section 2.2) are not shown, but both lie above 
the best-fit /3-model. 



for A2984 (compared to 1.3 ± 0.35; Table 3) and 2.6 ± 0.28 for 
AS84 (compared to 3.7 ± 0.96; Table 3). In both cases our masses 
are ~ 1 .4 times larger (albeit consistent to within the statistical un- 
certainties), which places an upper bound on uncertainties associ- 
ated with the approach employed. As a result, we conclude that our 
X-ray derived mass estimates (and corresponding rsoo and r25oo 
apertures) are reliable and that their measurement errors capture all 
reasonable sources of uncertainty. However, the impact of the mass 
measurement method is incorporated into our estimate of the sys- 
tematic error on gas fraction at the end of Section 5.4. The main rea- 
son for the systematic difference in our masses compared to those 
of Sivanandam et al. is a difference in the steepness of the outer 
gas density profile, which directly affects the pressure gradient and 
hence the enclosed hydrostatic mass, and we investigate this issue 
further in the next section. 



5.4 Choice of gas density model 

It is noteworthy that the issue of the gas density in the outer re- 
gions of galaxy groups and poor clusters is not well explored, but 
it is clearly of importance in measuring the gas mass. It is now 
well established that in more massive clusters the density profile 
beyond > 0.5 rsoo typically steepens relative to a power law ex- 
trapolated from smaller radii (i.e. a /?-model), from studies with 
Rosat (Vikhlinin, Forman, & Jones 1999), XMM-Newton (Croston 
et al. 2008) and Chandra (e.g. Ettori & Balestra 2009). However, no 
such steepening has been found in galaxy groups, where a /3-model 
provides a good description of the density profile out to rsoo (Ras- 



mussen & Ponman 2004). We therefore also fitted /3-models to the 
deprojected data points for all 5 clusters to cross-check the results 
from the AD08 cluster model, which does not use a /3-model for the 
gas density profile. The comparison is made in Fig. 1 3 and it is clear 
that for A2955, AS0296 & A2984 there is excellent agreement be- 
tween the two different models. Moreover, in each of these cases 
the data extend beyond rsoo, resulting in negligible differences in 
the gas mass within rsoo (< 3 per cent) between the two models. 
Beyond this radius the models tend to diverge, which reflects the 
lack of data points to anchor the fit, while differences between the 
two models inside of the innermost data point reflect the inability 
of the (flat cored) /3-model to capture the modest cuspiness in the 
density profiles. 

In making this comparison, it should be noted firstly that the 
ADOS model is jointly fitted to both /5gas(r) and T{r), assuming 
hydrostatic equilibrium and, also, that the ADOS model gas density 
profile is somewhat similar to a /3-model in the limit where there is 
no cool core component (see figure 1 of ADOS) - corresponding to 
t = 1, which is the case for A2955 (see Table 2 ). Furthermore, the 
ADOS model assumes a polytropic ICM, which may not represent 
the true physical state of the gas, although this assumption does 
not apply within the cool core (if one is present- i.e. t < 1 in 
Table 2). In the case of A2401 there is a difference between the 
extrapolated model profiles beyond the outermost data point, with 
the ADOS model steepening compared to the /3-model and yielding 
a lower gas mass within rsoo by 15 per cent. On the other hand, the 
two gas mass estimates agree to within 5 per cent for AS84, where 
extrapolation is also needed to reach rsoo. 
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It is clear from the 3 clusters with the best radial data cov- 
erage (A2984, A2955 & AS0296; Fig. 13) that no steepening of 
the gas density profile logarithmic slope is observed, within at least 
rsoo- This is quite different from the trend in more massive clusters 
and we discuss the implications of this behaviour in Section 5.7. 
However, our results are in accordance with Rasmussen & Pon- 
man (2004), who found that pgas(r') was well described by a /3- 
model over the entire radial range (which extended beyond ^rsoo) 
in two galaxy clusters of a similar mass scale (~2 kev). On the other 
hand. Sun et al. (2009) claim that the Vikhlinin et al. (2006) Pgas [r) 
model (which features a steepening logarithmic gradient) provides 
a good fit to their Chandra sample of 43 galaxy groups (although 
a number of these systems are not well covered azimuthally and/or 
radially beyond ^O.Srsoo)- 

Finally, as was the case with the total mass, we are also able to 
perform a cross-check of our gas mass results for A2985 and AS84, 
for which the same XMM-Newton data were analysed by Sivanan- 
dam et al. (2009), using a different analysis and modelling ap- 
proach. They generated XMM-Newton MOS X-ray surface bright- 
ness images, with corresponding exposure maps and background 
images and fitted these with a two-dimensional /3-model, convolved 
with the instrument point spread function. They then deprojected 
this fitted profile using the Abel transform, and applied a global 
spectral fit to determine the corresponding normalization in terms 
of gas density. The values of Mgas (< J"5oo) obtained by Sivanan- 
dam et al. are (in units of lO" Mg) 1.1 ± 0.08 for A2984 (com- 
pared to 1.21 ± 0.15; based on the data shown in Table 3) and 
2.3 ±0.13 for AS84 (compared to 2.77±0.61; Table 3), which are 
not significantly different from our own measurements. 

Although our estimates of the separate systematic effects of 
total and gas mass measurement method are not substantially dif- 
ferent from zero, they compound in the measurement of the gas 
fraction. Our best estimate of the systematic error in measuring 
/gas can be obtained from a comparison of our results with those 
of Sivanandam et al. (2009), who used different methods to us for 
determining both the total and gas mass (albeit applied to the same 
XMM-Newton data). The ratios of /gas within 7-500 measured with 
the AD08 model vs. Sivanandam et al. method are 0.83 ± 0.17 
and 0.85 ± 0.26 (for A2984 & AS84, respectively, with fractional 
errors summed in quadrature). While neither of these ratios is sig- 
nificantly different from 1, they are quite similar and somewhat less 
than one. Taking our ADOS model results as a baseline, this repre- 
sents a systematic increase of 15 per cent in /gas as a result of using 
a different analysis method. We therefore assign a systematic un- 
certainty of -1-15 per cent associated with the deprojection analysis 
and cluster modelling, but acting only in the direction of increas- 
ing /gas. This and the other contributions to the systematic error 
budget for the gas fraction are summarised in the middle column of 
Table 6, with negative values quantifying any potential decrease in 
/gas with respect to our baseline model. 

5.5 Other systematic effects 

The bootstrap re-sampling approach, used to evaluate uncertainties 
on all quantities derived from the cluster model, yields a direct es- 
timate of the statistical error that takes full account of intrinsic cor- 
relations (e.g between the gas and total mass in calculating /gas) 
and avoids the need for error propagation. However, it is important 
to consider sources of systematic uncertainty that could bias our 
results in a way that is not captured by these measurement errors. 

As already mentioned in Section 2.2, one source of systematic 
error is gas clumping (e.g. Urban et al. 2011; Nagai & Lau 2011, 





Percentage 


contribution 




f 

J gas 


J star 


Cluster modelling 


+15 /-O 




Non-hydrostatic equilibrium 


+0/-10 


+0/-4 


Outer bin density coiTection factor 


±7 




XMM-Newton calibration 


±5 




Luminosity function 




±8 


Mass to light ratio 




±15 


Total 


+27 / -22 


+23 / -28 



Table 6. Estimates of the positive and negative percentage contributions to 
the systematic errors on the gas and stellai' fractions. 

and references therein), which acts to bias the gas density higher. 
Since it is more prevalent in cluster outskirts, clumping would lead 
to a flatter density profile and a correspondingly overestimated gas 
mass and /gas, although this is difficult to quantify. In a similar 
fashion, non-thermal pressure support or merger-driven deviations 
from hydrostatic equilibrium could lead to an overestimated gas 
fraction, by ~5-20 per cent - the amount by which the total mass 
is likely underestimated (e.g. Kay et al. 2004; Rasia et al. 2006; 
Nagai et al. 2007; Piffaretti & Valdarnini 2008). In any case, both 
these phenomena would lead to problems in fitting the AD08 clus- 
ter model and likely result in a discrepant gas mass compared to 
that estimated from a surface brightness deprojection, using a /3- 
model for example. While we cannot rule this out, the indications 
from Section 5.4 are that any such differences are small compared 
to our measurement errors. 

To calculate the impact on /gas and /star of a +20 per cent 
systematic error on M500 resulting from non-hydrostatic bias, we 
adopt a similar approach to that detailed in Section 5.2. A new rsoo 
is calculated based on 1.2 x M500 (via Equation 2) and this aper- 
ture is used to measure the gas and stellar mass, which are then 
divided by 1.2 x Afsoo. In this case the stellar mass is calculated by 
simple extrapolation of the linear interpolation in log-log space of 
the spot values of /star measured at r25oo and rsoo. The resulting 
changes for the 5 clusters have mean values and standard errors of 
-(10 ± 0.6) per cent for /gas and -(4 ± 0.6) per cent for /star. Since 
this bias only acts one way (to cause the X-ray inferred mass to 
be underestimated), we take these mean values as the negative sys- 
tematic uncertainty on these quantities, with a corresponding zero 
contribution in the positive direction (see Table 6). 

Another source of systematic error is the model-dependent 
correction factor used to determine the gas density in the outermost 
bin (Section 2.2). In the case of A2955 & AS0296, we find that 
refitting the cluster model with this outermost bin omitted leads to 
a gas fraction within rsoo which is 6 per cent lower and 19 per 
cent higher, respectively. This is comparable with the fractional 
measurement errors (of 8 & 12 per cent) and, in view of the in- 
consistency of direction, we have assigned an average (symmetric) 
systematic error of ±7 per cent to account for this effect (i.e. the 
mean of -6 and +19, rounded up; see Table 6). 

Since our X-ray analysis is based entirely on XMM-Newton 
observations, we also need to consider the possible impact of sys- 
tematic calibration uncertainties. Although this is difficult to quan- 
tify, the most recent (Dec 2010) results of XMM-Newton cross- 
calibration with Chandra indicate*' that the effective area is subject 
to a ~5-10 per cent uncertainty, which would translate into a 2.5-5 



http : /xmm2 . esac . esa . int/external/xmm_sw_cal/ 
calib/ document at ion /index . shtml/ 
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per cent systematic uncertainty on the gas density/mass. However, 
this is based on older versions of the analysis software and the cal- 
ibration of both XMM-Newton and Chandra has improved since 
then. Nevertheless, we conservatively assume a (symmetric) sys- 
tematic error of ±5 per cent resulting from this (see Table 6). 

Sources of systematic error on the stellar fraction are some- 
what harder to quantify, but are dominated by the choice of mass- 
to-light ratio used to convert luminosity to stellar mass. As previ- 
ously mentioned, our chosen value of I band M/ L = 2.65 is some- 
what lower than the value of 3.6 used originally by GZZ07, but 
which now corrects for a 15 per cent dark matter contribution to 
the mass. The two main sources of systematic error in M/ L are the 
dark matter contribution and the faint end slope of the luminosity 
function, which we estimate to be ±15 and ±8 per cent, respec- 
tively (rightmost column of Table 6; Gonzalez et al, in preparation). 

Finally, it is worth noting that the impact of selection effects 
may be significant here, since the GZZ07 parent sample was se- 
lected on the basis of cluster optical properties and almost all large 
studies of cluster gas mass/fraction to date use X-ray selection. Se- 
lection effects are likely to lead to systematic differences in gas 
fraction, which will be biased higher for X-ray selected clusters; 
optically selected clusters could be biased low in /gas, if uncol- 
lapsed groups are included, although that is not the case for the 
five clusters considered here. On the basis of the weighted averages 
quoted in Table 5, the gas fraction for our (optically-selected) sam- 
ple is ~(19 ± 6) per cent lower than that predicted from the BCES 
regression to the X-ray selected literature sample at our correspond- 
ing mean M500. While optically-selected clusters (and particularly 
groups) often have at least a modest incidence of non-detections 
(e.g. Rasmussen et al. 2006; Balogh et al. 2011), it is interesting 
that the XMM-Newton follow-up of the GZZ07 sample has yielded 
good X-ray detections in all cases, with sufficient data quality to 
enable the gas fraction to be measured. However, we have ignored 
any systematic uncertainty associated with selection effects, since 
these are not relevant to the issue of baryon deficiency in individual 
clusters. 



5.6 Systematic error on the baryon fraction 

To assess the impact of systematic errors on the total baryon frac- 
tion within rsoo, we consider the separate contributions from the 
gas and stellar fractions. Based on our accounting of the sources of 
systematic uncertainty listed in Table 6, the most conservative ag- 
gregate values obtained (by a separate summation of the individual 
positive and negative contributions) are +211-22 per cent for /gas 
and -I-23/-28 per cent for /star - This assumes the most extreme case 
where contributions with opposite signs (e.g. the cluster modelling 
and non-hydrostatic biases) do not cancel each other out. 

The relative importance of both systematic and measurement 
errors on the measurement of the total baryon fraction can be 
gauged from Fig. 14, where the /star is plotted against /gas (both 
measured within rsoo). The measurement errors are shown as thin 
black lines, with systematic errors plotted in thick grey bars. The 
Universal mean /b and 'depleted cosmic' /b lie along diagonal 
lines in this plane, oriented at 45°, as a result of using an equal 
scale per unit length on each axis. Although the (conservative) sys- 
tematics dominate the measurement uncertainties for most of these 
systems, it is clear that at least 4 of the 5 clusters lie well below 
these limits. We can therefore be confident that a significant baryon 
deficit is present in these cases. 




Gas fraction within rsoo 

Figure 14. The location of the five clusters in the /star-/gas plane, com- 
pared to the Universal mean baryon fraction. The two different types of er- 
ror bar show the separate contributions from 1 a measurement (thin black) 
and conservative systematic (thick grey) uncertainties (see Table 6). 



5.7 Feedback and missing baryons 

There is considerable variation in /b for these five clusters, which 
results mainly from the scatter in /gas (Fig. 8, right panel), as the 
dominant baryonic component. While there are no indications of 
a systematic variation of /b with total mass (Fig. 9), the baryon 
fractions of these five clusters are lower than those of Gonzalez et 
al. (in preparation), which are generally more massive. Moreover, 
in the case of the 3 clusters where we can directly measure both 
gas and stars out to rsoo, there is evidently substantial baryon loss. 
However, it is not clear if this points to a trend with mass, or is the 
result of a break in self-similarity, as seen to occur at kT ~ 3.5 keV 
(corresponding to ~ 2.6 x lO^^M©) in a recent analysis of the X- 
ray luminosity-temperature relation (Maughan et al. 2012). 

It is well established that cosmological simulations without 
radiative cooling or galaxy feedback produce clusters that have 
baryon fractions within rsoo that are consistent with the 'depleted 
cosmic' /b level of ~0.9 (e.g. Ettori et al. 2006). The introduction 
of cooling allows baryons to transfer from the hot to the cold phase 
to form stars and leads to a systematic variation in gas fraction with 
halo mass (e.g. Muanwong et al. 2001). Indeed, Bode, Ostriker, & 
Vikhlinin (2009) point out that if stellar fraction varies as steeply as 
Afsoo"^ *^, then cooling alone can explain the trend of gas fraction 
with halo mass. However, radiative cooling produces no net loss 
of baryons- this is clearly inconsistent with our results, which also 
include the contribution from intracluster stars. Furthermore, the 
Bode et al. model assumes that the stars have the same radial distri- 
bution as the dark matter, whereas our results and those of GZZ07 
favour a less extended light distribution than that of the dark matter. 
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One possible explanation for the missing baryons is therefore 
galaxy feedback, for which the phenomenon of AGN outflows is 
the most likely candidate (e.g. Short & Thomas 2009; McCarthy 
et al. 2010, 201 1), whether prior to or after gravitational collapse. 
Indeed, it is now known that the feedback associated with galaxy 
formation and black hole growth is capable of altering the distribu- 
tion of both baryonic and dark matter even on large scales, thereby 
rendering cosmological inference sensitive to this localized bary- 
onic physics (Lueker et al. 2010; van Daalen et al. 2011). It is 
therefore of great importance to understand how cosmic feedback 
operates, and the corresponding balance of hot and cold baryons in 
groups and clusters of galaxies. Under this scenario, a diversity in 
feedback histories amongst the clusters could also account for the 
scatter in /gas, which dominates the variation in the total baryon 
content. 

Furthermore, the phenomenon of power-law-like pgas(j') in 
the outskirts of groups compared to a steepening in clusters is 
also seen in recent cosmological simulations where feedback has 
been incorporated and found able and necessary to produce galaxy 
groups that match many of the characteristics of observed sys- 
tems (McCarthy et al. 2010). For example, the shape of the median 
(group-dominated) gas entropy profile in figure 1 of McCarthy et al. 
(2010) roughly comprises a steeper inner power law (< 0.2r5oo), 
turning over to a flatter slope until ~0.9r5oo, before steepening 
again. This behaviour is a consequence of the increased impact of 
galaxy feedback in the shallower potential wells of groups com- 
pared to clusters which, firstly, displaces gas from inner to outer re- 
gions and, secondly, lowers the ICM pressure and therefore causes 
more gas accretion from outside the group (I. McCarthy, private 
communication). The net effect is a flattening of Pgas('") around 
rsoo in groups compared to clusters. 



6 CONCLUSIONS 

We have analysed the hot gas and stellar baryon content of five 
low-mass galaxy clusters drawn from the optically-selected sam- 
ple of GZZ07, which includes both a direct mapping of the in- 
tracluster light (ICL) as well as the galaxies' contribution, fully 
corrected for the effects of projection. We use hydrostatic X-ray 
cluster modelling of XMM-Newton observations to measure the gas 
total mass of these objects and the corresponding overdensity radii 
within which to perform an inventory of baryons (r25oo and rsoo). 
We summarize our main findings as follows. 

(i) For all five clusters our hydrostatic X-ray masses are well in 
excess of those originally estimated by GZZ07 on the basis of their 
galaxy velocity dispersions (a) and an assumed a — M500 relation, 
by a factor of between 1.7 and 5.7 times. We also calculate alter- 
native estimates of these X-ray masses using the M — Tx relation 
of Vikhlinin et al. (2006), and consider a range of possible expla- 
nations for this discrepancy (Sections 5.1 & 5.3). 

(ii) We study the stellar fraction vs. M500, compared to the other 
clusters with similar ICL mapping in GZZ07 and find that, despite 
the large change in total mass, the 5 clusters largely shift along the 
trend (since the higher mass implies larger rsoo, within which /star 
is measured). 

(iii) Within rsoo, we find weighted mean stellar (including ICL), 
gas and total baryon mass fractions of 0.026±0.003, 0.070±0.005 
and 0.096 ± 0.006, respectively, at a corresponding weighted mean 
M500 of (l.OSl^ ig) X M0 (see Table 5). This baryon fraction 
amounts to 57 per cent of the Universal mean, and implies a deficit 



of at least ~30 per cent compared to the 'depleted cosmic' /b level 
of ~0.9x fit /flm expected in the absence of cooling and feedback 
(e.g. Ettori et al. 2006). We consider the main sources of systematic 
error affecting our results and conclude that these are not sufficient 
to explain the observed shortfall in baryons in 4 of the 5 clusters. 

(iv) In the three cases where we can trace the X-ray emission out 
to rsoo, we find no evidence for steepening of the gas density pro- 
file, as commonly seen in more massive clusters. This behaviour in 
consistent with the increased effect of galaxy feedback in shallower 
potential wells seen in the OWLS cosmological simulations, which 
leads to gas accumulation around rsoo (see Section 5.7). This may 
also account for the baryon deficit observed within rsoo, and sug- 
gests that the missing baryons are located beyond this radius. 

As a result of the upward revision in mass, these systems are 
now no longer located as far towards the extreme low-mass end of 
the trend as previously thought (GZZ07), at a scale where the im- 
pact of cooling and feedback are most pronounced. Indeed, com- 
bined ICL photometry and X-ray mapping of many more low mass 
groups is clearly needed, to constrain the baryon balance in this im- 
portant regime. However, this remains observationally challenging, 
owing to the tension between the need for proximity (z <^ 0.1) in 
X-ray detection of faint group emission vs. the requirement of more 
compact angular size to provide stable and well-characterized back- 
grounds for diffuse optical light mapping, which is better achieved 
nearer z ~ O.L Given the uncertainly in modelling the gas density 
profile in the outskirts of groups and low-mass clusters, it is also 
very important to secure deep X-ray observations that can reach be- 
yond rsoo at high signal-to-noise ratio, to avoid extrapolation bias 
and thereby measure the hot gas mass directly. This would have the 
added benefit of probing a mass scale and spatial regions which are 
sensitive to cosmic feedback, and may help uncover at least some 
of the missing baryons. 
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